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1.0  Introduction 

This  report  presents  work  that  was  performed  under  the  Air  Force  Research  Laboratory’s 
Cooperative  Agreement  F30602-0 1-2-0540,  entitled,  “Simulation  of  Biomolecular 
Nanomechanical  Systems  (BioNEMS).”  Funding  for  this  effort  came  from  the  Defense 
Advanced  Research  Project  Agency’s  Simulation  of  Biological  System  (SIMBIOSYS)  Program. 
The  work  can  be  broadly  classified  into  three  tasks  around  which  the  report  is  organized.  In  the 
first  task,  the  goal  was  to  develop  computer  simulation  tools  and  theoretical  analyses  that  can 
guide  the  design  of  microdevices  that  sense  molecular  signatures  of  potential  pathogens  by 
translating  biomolecular  binding  to  mechanical  signals.  A  hierarchical  approach  was  adopted 
involving  close  synergy  between  experimental  and  computational  efforts,  and  models  were 
developed  for  phenomena  that  span  a  spectrum  of  time  and  length  scales  as  illustrated  in  Figure  1 
below.  The  results  of  this  work  were  then  used  to  establish  design  guidelines  and  utilized  in  the 
fabrication  of  a  prototype  device.  The  second  task  was  closely  related  to  the  first,  and  aimed  to 
develop  methods  to  detect  the  signals  generated  by  biomolecular  binding  using  electronics.  The 
third  task  explored  whether  new  discoveries  regarding  how  T  lymphocytes  of  the  immune 
system  detect  pathogens  can  be  exploited  to  create  synthetic  pathogen  detectors  that  exhibit 
extraordinary  sensitivity  and  selectivity.  This  work  examined  if  the  extraordinary  sensitivity  and 
selectivity  of  T  lymphocytes  for  antigen  could  be  maintained  if  the  antigen  presenting  cell  was 
replaced  with  a  synthetic  lipid  bilayer.  Computer  simulations  exploring  T  cell  signaling  were 
also  pursued.  The  work  performed  under  this  effort  will  be  presented  by  task. 
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Figure  1.  Illustration  of  simulation  organization  for  project. 
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2.0  Simulation  of  Biomolecular  Nanomechanical  Systems 


This  effort  was  among  the  first  to  show  that 
adsorption  of  biomolecules  on  one  surface  of  a 
microcantilever  generates  surface  stresses  that 
cause  the  cantilever  to  deflect,  and  the  adsorption 
of  a  second  molecule  that  binds  to  the  bound 
molecule  causes  the  deflection  to  change'’^.  This 
idea  defines  the  basic  platform  that  was  studied 
(Figure  2).  The  effect  has  been  demonstrated  for 
various  systems^.  The  goal  of  this  project  was: 

(i)  to  develop  simulation  methods  and  design 
rules  that  will  enable  the  design  of 
microcantilever-based  biosensing  systems;  and  (ii) 
to  develop  a  prototype  microdevice.  The  project 
involved  synergistic  theoretical,  computational  and 
experimental  work.  The  outcomes  of  both  the 
experimental  and  computational  efforts  are 
summarized  below. 

Theoretical  and  Computational  Effort 


♦  ^ 


Target  Molecule 


\ 


^Probe  Molecule 
•Gold 
SiN, 


/.♦/ 


Target  Binding 

^  Change  in 
deflection 


Figure  2.  Specific  biomolecular  interactions 
between  target  and  probe  molecules  alter 
intermolecular  interactions  within  a  self- 
assembled  monolayer  on  one  side  of  a  cantilever 
beam.  This  can  produce  a  sufficiently  large 
surface  stress  to  bend  the  cantilever  beam  and 
generate  motion. 


The  theoretical  and  computational  work  has  led  to  new  computational  tools.  More  importantly, 
the  theoretical  and  computational  efforts  identified  important  variables  that  affeet  mierodesign 
performanee  that  had  not  been  intuited  before  and  provided  design  rules  that  have  guided  the 
design  of  the  prototype  deviee  whose  commereial  applieations  are  currently  being  explored. 

This  theoretical  and  computational  work  led  to  aceomplishments  that  can  be  divided  into  four 
distinct  parts.  Each  of  these  is  deseribed  below. 

Nanomechanieal  forees  generated  by  surface  grafted  DNA 

Computational  and  theoretical  tools  were  first  developed  that  eould  caleulate  the  deflection  of 
microcantilevers  upon  adsorption  of  probe  molecules  and  changes  in  the  deflection  upon  binding 
of  a  target  biomolecule  to  the  layer  of  adsorbed  probe  moleeules.  In  short,  methods  were 
developed  that  allow  us  to  compute  how  binding  is  translated  into  a  mechanieal  signal,  and  to 
understand  how  various  conditions  and  the  identity  of  the  biomolecules  influences  eantilever 
defieetion.  The  focus  was  on  DNA  hybridization  as  a  prototype  system.  A  model  was 
developed  where  the  free  energy  eontributions  due  to  various  effects  (e.g.,  configurational 
entropy  of  adsorbed  molecules  and  counterions,  hydration  energies,  etc.)  were  earefully 
accounted  for.  Studies  revealed  three  important  points  that  had  not  been  understood  heretofore. 
First,  the  dominant  eontribution  to  eantilever  defieetion  is  hydration  forces,  not  configurational 
entropy  or  eleetrostaties.  Seeondly,  the  density  with  which  probe  molecules  are  grafted  on  to  the 
surface  strongly  influenees  cantilever  defieetion.  Thirdly,  the  cantilever  deflection  is  strongly 
dependent  upon  whether  or  not  the  grafting  of  probe  moleeules  on  the  surface  is  ordered  or  is 
disordered.  Disordered  distributions,  which  are  most  likely  in  praetice,  lead  to  much  larger 
cantilever  defleetions  (Figure  3).  These  findings  highlighted,  for  the  first  time,  that  molecular 
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level  self-assembly  controls  the  response  of  the  microdevice.  This  result,  and  associated  scaling 
relationships  were  crucial  for  the  experimental  program  (see  below)  which  ultimately  enabled  the 
development  of  a  prototype  device.  This  work  is  described  in  detail  in  an  article  published  in  the 
J.  Physical  Chemistry^. 


DNA  Grafting  Density,  a  (moiecuies/nm^) 


Figure  3.  Calculated  cantilever  deflection  as  a  function  of  the  grafting  density.  The  effects  of  disordered 
adsorption  highlighted,  for  the  first  time,  that  molecular  level  self-assembly  controls  microdevice  performance. 


Atomistic  understanding  of  DNA  hybridization  and  melting 

In  carrying  out  the  work  described  above,  it  was  discovered  that  potential  functions  describing 
the  hybridization  dynamics  of  DNA  were  not  well-developed.  The  reason  for  this  is  that  fully 
atomistic  (with  explicit  solvent)  simulations  of  DNA  hybridization  had  been  impossible  to  carry 
out  because  it  was  computationally  intractable  to  simulate  such  rare  events.  A  transition  path 
sampling  method^  was  applied  with  the  CHARMM  (Chemistry  at  HARvard  Macromolecular 
Mechanics/  program  for  macromolecular  simulations;  http://yuri.harvard.edu)  potential  to  enable 
the  study  of  DNA  hybridization.  This  first  of  its  kind  study  allowed  us  to  determine  the 
molecular  pathways  of  DNA  hybridization  and  melting  as  well  as  clarify  the  role  of  the  solvent 
in  these  processes.  Figure  4  shows  the  dominant  molecular  pathways  that  were  found.  This 
work  is  significant  for  two  other  reasons.  First,  it  provides  the  basic  information  necessary  for 
the  construction  of  accurate  potential  functions  for  future  use.  Secondly,  the  methodology  that 
was  developed  has  now  been  incorporated  in  the  CHARMM  program.  This  is  probably  the  most 
commonly  used  suite  of  programs  used  for  biomolecular  simulations,  and  so  our  contribution 
will  help  the  community  at  large.  This  work  was  been  published  in  the  Proceedings  of  the 
National  Academy  of  Sciences^. 
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Base  Rotates 


Figure  4.  Results  of  atomistic  molecular  dynamics  simulations  showing  the  major  routes  for  base-pair 
unbinding.  Details  regarding  the  free  energy  changes  are  described  in®. 


Understanding  the  dynamics  of  DNA  hybridization  on  cantilevers 


In  addition  to  the  thermodynamics  described  above,  cantilever  deflection  is  determined  by  the 
dynamics  of  DNA  hybridization  in  an  adsorbed  layer.  Adsorbing  the  probe  molecules  on  the 
surface  imposes  kinetic  constraints  on  the  hybridization  process  that  can  lead  to  incomplete 
binding  between  target  and  probe  molecules,  and  this  can  strongly  affect  the  response  of  the 
microdevice.  The  nucleation  of  DNA  hybridization  events  depends  upon  two  factors,  the 
probability  that  complementary  regions  overlap  and  the  rate  of  forming  a  nucleus  upon  overlap. 
A  stochastic  method  was  developed  to  investigate  how  the  structure  of  the  adsorbed  layer  of 
probe  molecules  on  the  cantilever  influences  the  transport  of  targets  on  to  the  layer  and 
subsequent  organization  of  complementary  regions  leading  to  nucleation  and  completion  of 
hybridization.  Specifically,  Master  Equations  that  can  be  solved  to  obtain  the  stochastic  events 
that  describe  probe  and  target  molecules  in  the  layer  were  formulated.  These  equations  were 
solved  to  obtain  the  initial  rate  of  hybridization.  For  the  grafting  densities  used  in  the 
experiments,  it  was  determined  that  transport  of  target  molecules  through  the  layer  is  not 
limiting.  Probe-probe  and  probe-target  interactions,  however,  decrease  the  number  of  available 
nucleation  sites.  This  effect  can  decrease  the  effective  hybridization  rates  by  over  an  order  of 
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magnitude  even  if  the  intrinsie  hybridization  kineties  is  unehanged  from  that  measured  in 
solution.  Higher  surfaee  eoverages  and  longer  probe  lengths  exaeerbate  these  effeets. 

Figure  5  shows  the  effeetive  rate  of  hybridization  as  a  function  of  the  grafting  density  of  probe 
molecules.  The  results  are  in  good  agreement  with  experimentally  measured  values  . 
Furthermore,  in  accordance  with  experiments^,  we  find  that  targets  that  preferentially  bind  to  the 
segments  near  the  grafted  end  of  the  probe  molecules  show  more  drastic  reductions  in 
hybridization  rates.  This  is  because  nucleation  sites  near  the  grafted  ends  are  least  accessible. 
These  results  suggest  that  grafting  density  must  be  controlled  during  assembly  of  the  probe 
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Experimental  data,  Henry  et  al.. 
Anal.  Biochem.  276,  204  (1999) 
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Figure  5.  The  predicted  dependence  of  the  effective  rate  of  hybridization  on  grafting  density,  and  comparison  with 

experiment. 


molecules  on  the  cantilever  if  reliable  microdevice  performance  is  to  be  obtained.  These  studies 
emphasized  the  importance  of  the  surface  chemistry  of  the  cantilever  for  designing  proper 
microdevices.  The  results  of  this  theoretical  work  have  been  reported  in  detail  in  a  publication  in 
the  J.  Chemical  Physics^. 

Molecular  Dynamics  Simulation  of  Hybridization  Kinetics  on  Cantilevers 


The  weakness  of  the  work  described  in  the  preceding  subsection  is  that  the  Master  Equations 
could  only  be  solved  easily  to  obtain  the  initial  hybridization  rates.  This  is  because  it  was  not 
possible  to  properly  account  for  changes  in  the  structure  of  the  probe  layer  during  hybridization 
using  this  framework.  So,  to  further  support  these  studies  emphasizing  the  importance  of 
grafting  density.  Molecular  Dynamics  simulations  were  carried  out.  These  simulations  were 
carried  out  using  standard  molecular  dynamics  with  the  potential  function  shown  in  Figure  6. 
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Molecular  dynamics  (MD) 

Potential:  modified  bead  &  spring  model 


V  =  v  +V  +  F 

non- bond  bond  angle 


i{r)  =  4f[(<T/r)'" 


Simulation  details 

NVT  ensemble;  T=300K;  PBC; 
Pg  =  0. 035-0. 35chains/nm2;one 
target;  each  probe/target:  20 
bps;  no  solvent;  10^  steps  ~  ps 

0.54nnn 
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Model  parameters* 

G  =  0.53  nm 
8  =  1 .67  kcal/mol 
m  =  250  grams/mol 
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F  =  42.3  kcal/mol 
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*Numerical  values  for  the  model  parameters  are  obtained  from  DNA  X-ray 
diffraction  and  bending  and  stretching  experiments  described  in  Ch.  3  and  5  of 
V.  A.  Bloomfiels,  D.  M.  Crothers,  and  I.  J.  Tinoco,  Nucleic  Acids:  Structures, 
Properties  and  Functions  (University  Science  Books,  Sausilito,  CA,  2000) 


Figure  6.  Potential  function  used  to  carry  out  Molecular  Dynamics  simulations. 


The  main  result  is  eneapsulated  in 
Figure  7  whieh  shows  how  at  larger 
grafting  density  hybridization  does  not 
occur,  but  rather,  the  target  chain  can 
bridge  across  more  than  one  probe 
molecule.  These  results  also 
suggested  scaling  relationships 
between  hybridization  efficiency  and 
various  parameters  that  control 
grafting  density  (e.g.,  salt 
concentration).  The  results  of  varying 
these  conditions  on  grafting  density 
and  cantilever  deflection  are  not 
detailed  here  because  they  are  reported 
in  the  description  of  the  experimental 
effort  (below)  as  the  computational  results  motivated  extensive  experiments  in  this  regard 
(reported  in  a  publication  in  Langmuir^). 

Experimental  Effort  on  Optical  Detection 

The  experimental  effort  in  this  project  was  originally  divided  into  two  components,  namely:  (i) 
development  of  a  microcantilever  array  chip  for  multiplexed  biomolecular  analysis;  (ii) 
measurement  of  surface  grafting  density  of  DNA.  Both  were  partially  motivated  by  results 
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Figure  7.  MD  simulations  for  typical  low  and  high  grafting 
densities.  For  high  grafting  densities,  bridging  is  a  severe 
problem. 
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emanating  from  our  computational  studies  of  microcantilever  deflection  at  equilibrium  upon 
DNA  hybridization.  In  addition  to  these  two,  another  component  to  the  experimental  work  was 
added,  namely:  (iii)  development  of  polymeric  membrane  based  biosensor.  Each  topic  will  be 
discussed  in  detail  now. 


Microcantilever  Array  Chip 

Figure  8  shows  a  section  of  the  5*  generation 
microcantilever  array  chip  that  was  developed''^.  This 
chip  contains  about  100  wells  with  each  well 
containing  4-5  cantilevers.  There  are  approximately 
500  cantilevers  on  each  chip.  Every  well  has  a  big 
I/O  for  injection  of  liquid  and  several  small  I/O’s  to 
prevent  bubbles  from  being  trapped  inside.  The 
design  of  the  fluidics  as  well  as  microcantilever 
structural  characteristics  was  done  in  close 
collaboration  with  CFDRC  researchers".  Each 
cantilever  beam  is  surrounded  by  its  own  reaction 
well  which  contains  two  fluidic  I/Os. 


Figure  8.  Optical  and  electron  micrographs  of 
5*  generation  cantilever  microarray  chip 
showing  the  whole  chip  and  a  single  well. 
Each  chip  contains  100  such  wells,  with  each 
well  containing  4-5  cantilevers.  All 
cantilevers  in  a  well  are  biofunctionalized  the 
same  way  through  fluid  inlef. 


The  microcantilevers  were  made  of  SiNx  and  fabricated  on  a  Si  wafer,  whereas  the  reaction  well 
is  defined  on  a  glass  wafer.  Each  cantilever  contains  a  large  paddle  for  reflecting  a  laser  beam 
that  is  used  for  the  optical  readout.  Figure  9  shows  the  schematic  diagram  of  the  optical  readout, 
where  a  collimated  beam  from  a  low-power  HeNe  laser  (1.5  mW)  was  expanded  and  reflected 
off  the  whole  microarray  chip  (see  Figure  9).  The  reflected  beam  from  the  whole  chip  was 
directed  towards  a  CCD  (charge-coupled  device)  camera.  Reflection  from  the  paddles  appear  as 
spots  on  the  CCD.  Because  of  the  initial  deflection  of  the  cantilevers,  the  reflections  from  the 
paddles  were  distinguished  from  the  reflection  off  the  silicon  surface.  Motion  of  these  CCD 
spots  was  monitored  to  measure  the 
deflections  of  all  the  cantilevers 
simultaneously.  Spots  from  500 
cantilevers  were  detected  and  their 
motion  tracked.  However,  lack  of 
automation  of  liquid  handling  and 
computation  power  have  forced  the 
experiments  to  be  conducted  using  40 
cantilevers  at  a  time. 


Figure  10  shows  plots  of  30  cantilever 
experiments  performed  almost 
simultaneously  (within  40  minutes)  in  8 
wells.  Figure  10  shows  cantilever 
deflection  as  a  function  of  time  for 
DNA  hybridization.  First,  non¬ 
complementary  target  ssDNA  20nt 
(nucleotide)  in  length  was  injected  into 


Actual  Image  Processed  Image 


Figure  9.  Schematic  diagram  of  the  optical  readout  system 
based  on  a  beam-deflection  technique  that  utilizes  a  single  laser 
and  a  CCD  camera  to  monitor  deflections  of  multiple 
cantilevers.  An  actual  CCD  image  of  reflection  spots  from  500 
cantilevers  is  shown.  Also  shown  is  the  processed  image, 
which  identifies  the  center  of  mass  (COM)  of  each  spot  and 
notes  the  COM  motion  for  each  spot  in  order  to  detect  the 
corresponding  cantilever  deflection. 
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all  the  wells  at  about  25  minutes.  This  resulted 
in  no  significant  cantilever  deflections. 
Complementary  target  ssDNA  20nt  in  length 
was  then  introduced  at  105  minutes  in  all  the 
wells.  The  cantilevers  immediately  responded 
and  produced  almost  similar  levels  of 
deflections.  In  all  cases,  the  deflections  were 
downward  away  from  the  gold  surface.  In  all 
cases,  the  cantilever  response  was  observed 
within  20-30  minutes,  which  is  much  shorter 
than  our  previous  observations. 

Achievements  and  challenges: 

1 .  Multiplexed  Experiments:  DNA 

hybridization  experiments  were  repeated 
successfully.  Drift  problems  were  resolved 
and  the  raw  data  usable,  as  shown  in  Figure 
10.  Intrawell  consistency  and  interwell 
repeatability  was  seen. 


0  30  60  90  120  150  180  0  30  60  90  120  150  180 


Time  Tmin^ 


2. 


Automation:  A  robotic  system  was 
designed  and  fabricated  and  used  to 
automatically  deliver  liquid  to  the 
microarray  in  a  precise  and  quick  manner. 


3. 


Figure  10.  Hybridization  Step:  Non-complementary 
DNA  was  injected  into  all  the  wells  around  25  minutes 
and  there  was  only  marginal  deflection. 
Complementary  DNA  was  injected  to  all  the  wells  at 
105  minutes  and  every  cantilever  started  to  deflect 
away  from  the  gold  surface.  Spikes  between  the  non- 
complementary  and  complementary  DNA  injection 
were  buffer  injections. 


Surface  Chemistry:  A  signifieant  amount 
of  time  and  effort  was  devoted  to 
developing  and  integrating  surfaee 

chemistry  to  bind  receptor  molecules  and  prevent  non-speeific  binding.  The  options  tried  for 
binding  proteins  to  the  gold  coated  cantilever  were:  (i)  DTSSP-NHS 
(Dithiobis(sulfosuccinimidylpropionate)  -  Nance-Horan  syndrome)  /  protein  using  amine 
chemistry;  (ii)  thiol-biotin/streptavidin/biotinylated  protein;  and  (iii)  thiol-alkane- 
NHS/protein  using  amine  ehemistry.  To  reduce  non-specifie  binding,  polyethylene  glyeol 
(PEG)  was  attached  to  the  bare  nitride  cantilever  surfaee  using  silane  chemistry.  The  pros 
and  eons  of  these  options  were  evaluated  and  eurrently  options  (i)  and  (iii)  are  generally 
used. 


4.  The  cantilever-based  antibody  array  was  used  for  three  types  of  experiments:  (i)  quantitative 
analysis  of  proteins;  (ii)  direet  deteetion  of  pathogens  sueh  as  E.  coli  by  antigen-antibody 
binding  to  a  membrane  protein;  and  (iii)  quantitative  analysis  of  enzyme  activity. 
Preliminary  results  in  both  experiments,  ineluding  controls,  were  positive. 


8 


Surface  Grafting  Density  of  Probe  Molecules 


Measurement  Technique  and  DNA  Measurements 

The  effort  to  measure  ssDNA  surface  grafting  density  was  in  direct  response  to  the  theoretical 
effort  whieh  showed  that  eantilever  defleetions  depend  exponentially  on  grafting  density.  This 
was  not  known  before  the  beginning  of  this  projeet.  The  teehnique  shown  in  Figure  1 1  was  used. 
Fluorescently-labeled  single-stranded  DNA  (ssDNA)  were  first  grafted  onto  a  Au-eoated  surfaee. 
They  are  then  reacted  with  alkane  thiols,  whieh  eompetitively  bind  with  the  Au  surfaee  and 
remove  the  ssDNA  into  solution.  The  concentration  of  ssDNA  in  solution  was  measured  either 
through  UV  absorption  or  through  fluorescenee.  By  measuring  the  ssDNA  coneentration  in 
solution,  fluorescence  quenehing  by  gold  and  fluoreseenee  eoupling  between  neighboring 
ssDNA  were  eliminated. 

A  mercury  lamp  was  used  for  exeitation  and  a  slow-sean  CCD  camera  for  detection  of 
fluorescence.  The  experiments  for  different  DNA  lengths  and  salt  concentrations  were  repeated 

Fluorescent 


Surface  Grafted  Thiol  Fluorescence/ 

single-stranded  DNA  Absorption 

Measurement 


Figure  11.  Method  to  measure  surfaee  grafting  density  of  ssDNA  onto  a  gold  surface. 

and  the  results  are  shown  in  Figure  12.  These  results  eorrespond  well  to  the  results  obtained 
earlier  using  the  laser-photodiode  approach.  Previous  experiments  indicated  that  the  interesting 
range  of  salt  coneentrations  in  whieh  the  grafting  density  ehanges  appreciably  lies  from  10-100 
mM  and  therefore  the  research  focussed  on  this  range  of  concentrations.  The  results  of  this  work 
was  published  in  a  paper  that  appeared  in  Langmuir'^. 

Grafting  Density  of  Proteins  on  Gold  Chips 

An  effort  in  measuring  the  grafting  density  of  antibodies  on  a  surface  was  completed'^.  The 
surface  density  of  anti-PSA  (prostate-speeifie  antigen)  on  gold  was  determined  using  the  same 
fluoreseenee  teehniques  deseribed  earlier.  The  protein  was  attached  to  the  gold  surface  using  a 
hetero-biflmetional  erosslinker  (DTSSP).  The  immobilized  protein  was  then  displaeed  into 
solution  by  adding  mercapto-ethanol,  which  displaces  the  DTSSP.  Finally,  the  displaced  protein 
in  solution  was  quantitated  using  the  CCD  eamera  and  calibrated  against  a  set  of  protein 
solutions  of  known  eoncentrations.  The  surfaee  density  of  the  immobilized  antibody  was 
eomputed  as  0.0055  nm'^,  whieh  roughly  corresponds  to  about  25%  surface  coverage  assuming  a 
typieal  size  of  an  IgG  (immunoglobin)  molecule. 
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Figure  12.  Grafting  density  of  ssDNA  measured  as  a  function  of  ssDNA  length  and  salt  concentration. 


DNA  Hybridization  Efficiency 

Based  on  the  discussions  at  a  SIMBIOSYS  Principal  Investigator’s  Meeting  (Sept.  2003  in 
Monterrey,  CA),  experiments  were  performed  using  the  above-mentioned  approach  to  study 
DNA  hybridization  efficiency  as  a  function  of  DNA  length,  salt  concentration  and  base-pair 
mismatches.  The  data  in  Figure  13  shows  that  efficiency  decreases  dramatically  with  increasing 
salt  concentration,  presumably  because  of  molecular  crowding  (increasing  grafting  density  - 


Figure  13.  DNA  hybridization  density  and  efficiency  as  a  function  of  DNA  length  and  concentration.  The  DNA 
strands  are  fully  complemenfary  in  these  experiments. 


Figure  12).  However,  increasing  grafting  density  and  decreasing  hybridization  efficiency  almost 
cancel  each  other  producing  an  almost  constant  hybridization  density.  Figure  14  shows  the  data 
for  hybridization  density  as  a  function  of  base  pair  mismatches.  The  data  suggests  that 
mismatches  at  the  ends  do  not  alter  density,  while  those  in  the  middle  have  a  marked  effect. 
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Figure  14.  DNA  hybridization  density  as  function  of  salt  concentration  for  fully  complemenfary,  single  base  pair 
mismatch  at  the  proximal  3’  end,  single  base  pair  mismatch  in  the  middle,  and  triple  base  pair  mismatches  in  the  middle. 


Antigen-Antibody  Binding  Efficiency 

Experiments  were  performed  to  detect  the  degree  of  antigen-antibody  binding  using 
fluorescence.  The  purpose  of  the  experiment  was  to  test  the  effect  of  antibody  immobilization. 
Two  approaches  were  used:  (i)  NHS-SS-Biotin  (biotin  disulphide  N-hydroxysuccinimide  ester) 
/Streptavidin/Biotinylated-AntiPSA  bound  to  PSA;  (ii)  DTSSP-AntiPSA  bound  to  PSA.  For  the 
NHS-SS-Biotin  approach,  the  bound  PSA  density  was  found  to  be  0.0046  nm'^  whereas  for  the 
DTSSP  approach,  the  density  was  found  to  be  0.0065  nm'^.  Since  the  typical  size  of  IgG  is  about 
5  nm,  the  binding  density  suggests  a  surface  coverage  of  20-25  %. 

The  influence  of  grafting  density  on  the  cantilever  deflection  signal  was  seen  by  using  the 
control  in  grafting  density  of  DNA  molecules^.  Figure  15  shows  that  a  scaling  relation  can  be 
obtained  for  cantilever  deflection  by  plotting  it  as  a  function  of  hybridization  density. 
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Figure  15.  Grafting  density  of  probe  DNA  and  a  hybridization  density  of  target  DNA  as  a  flinetion  of  salt 
eoneentration.  This  can  be  plotted  as  hybridization  efficiency  as  a  function  of  grafting  density  and  finally  cantilever 
deflection  as  a  function  of  hybridization  density. 


Using  the  surface  chemistry  that  we 
have  developed,  we  were  able  to  assay 
at  different  concentrations.  Figure  16 
shows  our  measurement  of  PSA  at 
concentrations  down  to  1  ng/ml,  against 
a  background  of  BSA  (bovine  serum 
albumin)  at  2  mg/ml We  are 
currently  trying  to  detect  such  proteins 
in  complex  mixtures  such  as  serum. 

Enzyme  Activity 


The  activity  of  a  kinase  during 
phosphorylation  of  a  serine  residue  of  a 
peptide  attached  to  a  cantilever  was 


fPSA  concentration  (ng/ml) 

Figure  16.  Measurement  of  PSA  against  a  background 
of  BSA  at  2  mg/ml. 


detected*^.  Figure  17  shows  the  schematic. 

A  monolayer  of  a  peptide  from  p53  protein  containing  a  recognition  element  for  the  DNA 
protein  kinase  and  serine  phosphorylation  site  was  attached  to  the  cantilever  using  thiol 
chemistry.  Then  a  DNA  protein  kinase,  along  with  its  reagents,  were  introduced  in  the  reaction 
well.  The  data  in  Figure  17  suggests  that  when  ATP  (adenosine  triphosphate)  is  introduced,  the 
cantilever  bends  significantly,  indicating  phosphorylation.  This  was  verified  by  the  experiment 
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with  ATP-analog,  which  cannot  be  hydrolyzed,  which  shows  no  signal.  In  addition,  when  a 
kinase  inhibitor  was  used,  there  was  no  signal  as  well.  This  clearly  proves  that  phosphorylation 


Figure  17.  Schematic  showing  p53  peptide  on  the  cantilever,  and  subsequent  phosphorylation  by  DNA  protein 
kinase.  The  cantilever  deflection  and  reaction  rates  are  distinctly  different  when  ATP  was  used  as  opposed  to 
when  ATP  analog,  or  an  inhibitor  was  used. 

can  be  studied  in  real  time  to  obtain  kinetics. 

Polymer-Membrane  Based  Biosensors 


A  new  experimental  effort  that  was  not  part  of  the  original  proposal  was  developed  for  several 
reasons: 

(A)  The  elasticity  modulus  of  polymers  is  two  orders  of  magnitude  lower  than  most 
microfabrication  materials  such  as  SiNx,  Si,  etc.  This  makes  the  mechanical  device  more 
compliant. 

(B)  Polymers  like  parylene  are  bio-compatible  and  inert. 

(C)  Polymers  are  hydrophobic,  which  can  allow  device  design  that  enables  selecting 
functionalization. 

(D)  Easier  device  fabrication  may  be  possible. 

(E)  Membranes  are  fixed  at  both  ends  and 
hence,  may  have  less  drift. 

A  novel  parylene  micro-membrane  surface 
stress  sensor  array  using  capacitive  signal 
was  developed  (see  Figure  18).  This  sensor 
exploits  the  low  elasticity  modulus,  better 
chemical  resistance  and  biocompatibility  of 
polymers  like  parylene  when  compared  to  the 
silicon  based  materials  used  in  traditional 
microfabrication.  The  salient  features  of  the 
sensor  are  that  it:  (i)  is  label  free;  (ii)  is  a 
universal  platform  -  suitable  for  both 
chemical  and  biological  sensing;  (iii)  uses 
electronic  (capacitive  detection)  readout;  (iv) 
has  integrated  microfluidics  for  addressing 
individual  sensors  on  the  chip;  (v)  is  capable 


1 - 

_ 

Figure  18.  Binding  of  analyte  species  to  the  probes 
immobilized  on  the  membrane  surface  changes  the 
surface  stress.  The  compressive  surface  stress  change 
as  shown  in  this  figure  causes  an  out  of  plane 
deflection  that  increases  the  air-gap  and  reduces  the 
capacitance. 
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of  handling  both  liquid  and  gas  samples;  (vi)  is  made  using  standard  low  temperature 
microfabrication  processes  (<  120  °C);  and  (vii)  can  readily  be  scaled  and  multiplexed. 

Finite  element  modeling  was  used  for 
optimization  of  sensor  geometry 
parameters.  A  microfabrication 
process  sequence  that  is  compatible 
with  parylene  processing  was 
designed  using  surface  and  bulk 
microfabrication  techniques  to 
manufacture  the  prototype  of  the 
sensor  array.  Figure  19  shows  optical 
micrographs  of  the  parylene  surface 
stress  sensors  with  integrated 
capacitance  measurements.  This 
sensor  array  platform  can  be  tuned  for 
sensing  many  different  targets  by 
using  suitable  sensor  coatings  and 
these  coatings  play  an  important  role 
in  determining  the  sensitivity  and 
resolution  of  the  sensor.  In  the  current 
work,  alkanethiols  with  different  functional  end  groups  were  chosen  as  sensor  coating  materials. 
Alkanethiol  molecules  have  a  strong  affinity  for  gold  surfaces  and  this  property  was  exploited  for 
coating  formation  on  the  gold  electrodes  of  the  membrane  sensor.  Furthermore,  the  use  of 
alkanethiol  coating  required  minimal  modification  in  the  sensor  microfabrication  process 
sequence. 

Figure  20  shows  the  results  of  multiple  functionalized  membrane  sensors  exposed  to  toluene  and 
isopropyl  alcohol.  It  is  clear  that  the  different  surface  coatings  produce  different  signals. 
However,  the  simple  alkane  thiol  molecules  do  not  have  selectivity  to  distinguish  between 
mixtures  of  gases. 


Figure  19.  (a)  The  micro-fabricated  sensor  chip  in  comparison 
with  a  nickel;  chip  has  a  3X3  array  of  sensors  unit  cells;  (b)  chip 
with  bonded  microfluidic  cover;  (c)  unit  cell  with  layout  for 
single  ended  capacitance  measurement;  (d)  unit  cell  with  layout 
for  differential  capacitance  measurement;  (e)  close-up  view  of  the 
released  membrane  sensor. 


Time  (seconds)  Time  (seconds) 

Figure  20.  Normalized  capacitance  vs.  time  on  (a)  toluene  and  (b)  isopropyl  alcohol  exposure  for  sensor 
membranes  coated  with  different  thiol  molecules.  The  terminal  groups  were:  MUA  -COOH;  MUO  -  OH;  DOT 
-CH3. 

These  results  now  appear  in  a  paper^"^  and  a  PhD  thesis^^. 
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The  activity  under  this  portion  of  the  research  accomplished  the  following: 

1.  Creation  of  a  new  electronic  method  of  detecting  surface  stress  using  parylene 
membranes  and  capacitance  detection. 

2.  Development  of  a  new  fabrication  process. 

3.  Concept  fabrication  and  performance  testing. 

4.  Demonstrated  chemical  sensing  using  parylene  membrane  devices. 

5.  Establishment  of  a  room  temperature  adhesive  bonding  process  which: 

a.  Designed  an  experiment  for  determining  the  optimum  process  parameters. 
Substrates  tested:  Silicon,  Glass,  PDMS,  Parylene 

Adhesives:  UV  curable  adhesive,  uncured  PDMS 
Burst  pressure  chosen  as  device  failure  criteria. 

b.  Completed  the  test  fixture  and  test  device  fabrication. 

c.  Completed  a  portion  of  device  testing.  Microfluidic  channels,  made  by  bonding 
glass  and  silicon  wafers  using  UV  curable  adhesives,  could  withstand  200  psi 
pressure  without  leaking. 

Figure  21  shows  the  bonding  process  while  Figure  22  shows  the  results  of  burst  pressure  as  a 
function  of  size  of  a  square  cavity*^.  It  is  evident  that  the  bonding  process  can  withstand 
pressure  beyond  100  psi,  which  is  adequate  for  many  biofluidic  chip  applications. 


r 


(b)  (d) 


Figure  21.  Adhesive  transfer  process;  (a)  uniform  adhesive  layer  on  a  wafer  obtained  by  spinning,  (b)  the  cover  wafer  with  patterns 
brought  into  contact  with  the  adhesive,  (c)  selective  transfer  of  the  adhesive  to  the  cover  wafer,  (d)  the  cover  wafer  aligned  with  a 
device  wafer,  brought  into  contact  and  bonded. 
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Figure  22.  The  burst  pressure  variation  with  square  device.  Each  point  represents  one  device.  Arrows 
indicate  that  the  bonding  did  not  fail  at  the  indicated  pressure. 
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Desired  performance  metrics  for  microdevice 


Label-free  detection 

No  reporter  molecules;  Reduced  reagent  cost 
Simultaneous  detection  of  multiple  pathogens,  molecules,  or  substances 
Fast  (<  5  min)  detection  time  after  sample  preparation 
Versatile  platform  for  detection 

Nucleic  acids;  antigen-antibody;  gas  or  liquid-based  chemicals/toxins 
Built-in  calibration  for  specific  environments 

High  specificity  and  sensitivity 

Presence  or  absence  of  75bp  (base  pairs)  DNA  sequence 

Ideally  lOfg  of  DNA  from  viruses  and  lOOfg  (femtogram)  of  DNA  from  bacteria 

Hand-held  battery-operated  device 
Commercialization  and  Technology  Transfer 


Based  on  the  technology  developed  under  this  program,  a  new  company,  Kalinex  Inc,  was 
launched  to  commercialize  this  technology  for  chemical  and  biological  sensing.  Kalinex 
received  a  small  seedling  from  DARPA  and  has  also  raised  private  investment  from  individuals 
and  venture  capital.  Kalinex  currently  employs  three  full  time  employees  and  three  consultants, 
and  is  in  the  process  of  making  a  commercial  product  for  homeland  security  and  industrial 
applications.  Kalinex  has  also  received  the  option  to  exclusively  license  the  patents  that  were 
filed  by  UC  Berkeley  (see  Appendix)  based  on  the  research  conducted  under  this  program. 
Kalinex  is  holding  discussions  with  private  companies  as  well  as  with  Transportation  Security 
Agency  for  launching  this  product. 


3.0  Electronic  Detection  on  MOSFET-Embedded  Microcantilevers 

The  objective  of  this  task  was  to  develop  an  electronic  detection  concept  for  integrated  remotely- 
addressable  bio-chem  sensors,  based  on  MOSFET-embedded  microcantilevers.  The  task 
successfully  demonstrated  the  core  objective  and  validated  with  peer-reviewed  archival 
publications.  The  research  also  resulted  in  intellectual  property  (IP)  related  to  MOSFET  sensing 
approach,  and  a  patent  application  has  been  filed  on  behalf  of  Northwestern  University  (NU). 

The  specific  project  goals  and  research  objectives  include: 

1)  Simulation  of  microcantilevers  to  optimize  high  surface  stress  region, 

2)  Design  and  layout  of  MOSFET-embedded  microcantilevers  for  biomolecular  detection, 

3)  Fabrication  of  prototype  MOSFET  on  cantilever,  and  preliminary  testing  for  performance 
evaluation 

The  project  approach  is  schematically  presented  in  Figure  23. 
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Figure  23.  Schematic  illustration  of  MOSFET  cantilevers  approach. 

Simulation  of  Microcantilevers  for  Optimization 


To  properly  identify  the  high  surfaee 
stress  region  on  the  eantilever,  the 
meehanical  properties  were  modeled 
and  simulated  using  ANSYS®  and 
Intellisuite®  finite-element  analysis 
software.  Both  packages  are  available 
at  NU  and  have  been  used  by  a  team 
of  experienced 

microelectromechanical  systems 


(MEMS)  researchers  for  cantilever 
devices.  Solid  models  of  the 
cantilevers  were  constructed  and 
proper  mesh  strategies  were 

determined  to  simulate  cantilever  performance  (e.g.,  Figure  24). 


Figure  24.  Top  view  of  a  Si  multilayer  cantilever  with  a  20pm 
anchor  area  showing  stress  distribution. 


Optimization  of  cantilever  width,  thickness  and  resonance  frequencies  was  carried  out  with  the 
software  packages  mentioned  above.  Specific  cantilever  designs  were  targeted  to  optimize  stress 
localization  at  the  base  of  the  cantilever  for  subsequent  measurement  via  MOSFET-stress 
transistors. 
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In  addition,  a  resonance  frequency  analysis  of  various  cantilever  designs  was  carried  out  to 
analyze  and  optimize  effects  of  mass  loading.  The  FEA  simulations  reduce  the  need  to  perform 
detailed  experimental  optimization  of  the  device  characteristics.  Numerous  simulation  using 
different  cantilever  geometries  and  specifications  were  carried  out.  The  lengths  of  the  cantilevers 
in  our  simulations  ranged  from  200-450  pm  and  the  thickness  from  1. 5-3.0  pm.  The  maximum 
stress  was  found  to  be  generated  at  the  base  of  the  cantilever.  This  area  is  referred  to  as  the  stress 
concentration  region  (SCR).  These  simulated  results  helped  us  in  guiding  where  to  locate 
MOSFETs  on  the  cantilever. 

MOSFET  Device  Fabrication  and  Testing 


Layout  of  the  on-chip  MOSFET  transistor  on  the  cantilever  was  carried  out  using  L-EDIT  design 
CAD.  Seven  masks  were  designed  using  L-EDIT  to  accomplish  this  goal  and  mask  layout  using 
L-EDIT  was  converted  into  a  pattern  generator  file  using  the  PG  generation  system  at  the 
University  of  Illinois  at  Chicago  (UIC)  microfabrication  facility.  Seven  masks  on  quartz  glass 
plates  were  developed.  The  liquid  cell  around  the  cantilever  chip  was  designed  and  laid  out  for 
testing  in  liquid  environment. 

Si  cantilever  fabrication  was  carried  out  on  a  SOI  (silicon  on  insulator)  wafer  with  a  device  layer 
of  10  pm.  Standard  Si  technology  processes  such  as  lithography,  PECVD  (plasma-enhanced 
chemical  vapor  deposition),  metal  deposition,  etching,  oxidation  and  diffusion  process  were 
carried  out  to  fabricate  the  micro  cantilever  with  embedded  MOSFET  on  the  high  stress  region 
of  the  cantilever.  Most  of  the  microfabrication  was  carried  out  in  the  microelectronics  laboratory 
at  UIC.  Ion-implantation  for  source  and  drain  were  carried  out  at  the  CNF  (Cornell  NanoScale 
Science  &  Technology  Facility)  user  facility.  Transistors  were  placed  along  the  longitudinal 
direction  of  the  cantilever  where  stress  is  maximum  and  a  differential  approach  was  used  to 
measure  current  change  as  a  result  of  cantilever  bending. 

Scanning  electron  microscopy  (SEM)  image 
of  one  pair  of  identical  cantilevers  from  an 
initial  50  x  1  microcantilever  array  is  shown 
in  Figure  25.  The  pair  consists  of  one 
microcantilever  coated  with  a  thin  film  of 
Cr/Au  for  immobilization  of  probe 
molecules,  typically  with  thiol  chemistry, 
and  the  other  was  uncoated  and  acted  as  the 
reference.  The  differential  drain  current 
between  the  sensing  and  the  reference 
microcantilevers,  which  further  minimizes 
systematic  noise  and  environmental 
perturbations,  forms  the  basis  for  the 
MOSFET  electronic  detection  (see  Figure 
23).  The  differential  signal  can  be  fed  into  a 
CMOS  based  differential  amplifier  for 
electronic  readout  at  the  chip  level  for  future 


Figure  25.  SEM  image  of  pair  of  MOSFET 
cantilevers,  indicating  contacts. 
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development.  MOSFET-embedded  micro-cantilevers  with  the  thickness  of  1.5  to  2pm  and  length 
ranging  from  200  to  300pm  were  fabricated,  with  the  separation  of  about  250  pm  between  the 
reference  and  sensing  microcantilevers  for  simplicity.  The  transistor  was  located  about  2  to  4  pm 
from  the  cantilever  base  and  the  W/L  ratio  of  source  and  drain  was  around  10  to  achieve  high 
transconductance.  The  resonance  frequency  of  the  MOSFET-embedded  microcantilevers  was 
around  100  to  150  kHz.  Each  array  was  designed  to  have  identical  sensor  (i.e.,  Cr/Au-coated) 
and  SiNx  reference  microcantilevers  for  differential  output  to  minimize  systematic  noise  and 
possible  false  positives.  The  residual  stress  that  may  be  introduced  by  applying  a  thin  layer  (30 
nm)  of  gold  coating  on  one  side  of  the  cantilever,  does  not  create  any  notable  difference  in 
MOSFET  current-voltage  characteristics  when  compared  with  those  of  Si3N4  reference 
cantilevers. 

The  drain  current  sensitivity  of  MOSFET-embedded  microcantilevers  to  bending  was  validated 
using  a  high-resolution  nanomanipulator  with  <5  nm  vertical  (Z-axis)  resolution.  A  decrease  in 
drain  current  between  0.1  to  0.2  mA  per  nanometer  microcantilever  deflection  was  observed. 
Changes  in  the  drain  current  by  almost  one  order  of  magnitude  between  few  to  150  nanometer  of 
microcantilever  bending  were  noted.  The  bending  results  indicated  that  the  MOSFET  deflection 
sensitivity  was  of  the  same  order  as  optical  detection  as  inferred  from  the  prior  literature  and 
almost  one-to-two  order  of  magnitude  higher  than  existing  active  and  passive  detection 
technologies. 


3 
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Figure  26.  Noise  and  S/N  data  from  MOSFET  cantilevers. 


Moreover,  MOSFETs  have  large  signal-to-noise  (S/N)  ratios  caused  by  the  large  change  in  drain 
current  relative  to  the  concomitant  small  noise  density  in  Figure  26.  The  three  curves  show  1/f  or 
flicker  noise,  the  dominant  source  of  noise  in  MOSFETs  at  low  frequencies,  at  three  different 
voltages  (23,24).  The  current  noise  of  40-60  nA,  which  was  calculated  by  integrating  the  spectral 
power  density  over  1/f  bandwidth  for  different  gate  voltages,  was  significantly  lower  than  that 
reported  for  piezoresistive  based  microcantilever  based  sensors  (2.7  pV)  (25).  Given  that 
MOSFET  current  sensitivity  was  around  0.1  to  0.2  mA/nm  of  cantilever  deflection,  low  detection 
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limits  can  be  readily  achieved  with  large  S/N  ratio.  The  noise  density  could  likely  be  further 
reduced  in  subsequent  generation  of  these  devices  by  standard  processing  steps,  involving 
optimization  of  doping  concentration  and  minimizing  the  interface  traps. 

Surface  Immobilization  of  Receptors 

For  biomolecular  binding  experiments,  the  MOSFET-embedded  microcantilevers  were  cleaned 
sequentially  in  acetone,  isopropanol-2,  and  methanol  for  10  min  each,  followed  by  UV-cleaning 
for  25  min,  and  functionalized  with  DTSSP  ([3,3'-Dithiobis(sulfosuccinimidylpropionate)], 
Pierce  Chemical  Company),  a  linker  molecule  involved  in  immobilizing  streptavidin  and 
antibodies  to  the  gold-coated  microcantilever  surface.  Streptavidin  (Pierce)  was  subsequently 
immobilized  on  the  microcantilever  surface  by  incubating  overnight  in  a  10  qg/ml  streptavidin 
solution  prepared  with  phosphate  buffered  saline  (PBS,  pH  =  7.4).  This  immobilization  method 
provides  a  tight  streptavidin  layer  with  uniform  density  on  gold  for  efficient  binding  of  biotin. 
All  of  the  non-specific  binding  sites  were  blocked  with  bovine  serum  albumin  (BSA).  For 
detection  experiments,  the  functionalized  microcantilevers  were  exposed  to  100  fg/ml,  100 
pg/ml,  and  100  ng/ml  of  target  biotin  in  PBS. 


Electronic  Signature  of  Microcantilever  Binding 


MOSFET  transistors  were 
passivated  with  silicon 
nitride  thin  coating  (30  nm) 
and  electrical  contacts  were 
isolated  for  the  binding 
measurements  in  the  fluidic 
environment.  The  measured 
Id  versus  Vds  characteristics 
for  n-MOSFET-embedded 
transistor,  at  Vg=5V,  shows 
a  negligible  change  in  Id 
(Figure  27 A)  when  the 
streptavidin  immobilized 
gold  microcantilevers  are 
immersed  in  PBS. 

Microcantilever  bending  as 
a  result  of  streptavidin- 
biotin  binding  leads  to 
decreases  in  Id  as  the 
concentration  of  biotin 
increases  from  100  fg/ml  to 
100  ng/ml.  No  drain 

current  change  was  seen  in  SiNx  cantilevers  with  biotin,  where  no  binding  events  occurred. 
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Figure  27A.  MOSFET  characteristics  for  Streptavidin-Biotin  binding 
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Similar  experiments  were 
performed  for  deteetion  of 
goat  antibodies  (secondary 
IgG)  by  rabbit  antibodies 
(primary  immunoglobulin  G, 

IgG)  with  the  embedded 
MOSFET.  After  the  cleaning 
procedure,  the  MOSFET- 
embedded  microcantilevers 
were  first  functionalized  with 
DTSSP  as  a  linker,  and 
incubated  overnight  in  0. 1 
mg/ml  rabbit  anti-goat  IgG 
(Pierce  Chemical  Company) 
prepared  in  PBS  for 
immobilization.  BSA  was 
again  used  as  an  agent  to 
block  non-specific  binding 
sites.  The  functionalized 
microcantilevers  were 
exposed  to  O.lmg/ml  of  goat 
anti-rabbit  IgG  in  PBS  for 
binding  experiments. 

The  measured  Id  versus  Vds  characteristics  Vg=  5  V  (Figure  27B)  again  shows  no  change  in  the 
drain  current  when  SiNx  and  rabbit  IgG  coated  cantilevers  are  immersed  in  PBS.  When  O.I 
mg/ml  of  goat  anti-rabbit  IgG  was  introduced,  an  almost  two  orders  of  magnitude  change  in  Id 
was  observed,  indicative  of  microcantilever  bending  as  a  result  of  antibody-secondary  antibody 
binding.  The  SiNx  reference  cantilever  remained  the  same  after  injecting  the  target. 

Thus,  a  novel  label-  and  optics-free  detection  platform  for  biomolecular  binding  was  developed. 
It  involves  high  sensitivity  detection  of  cantilever  bending  via  embedded  MOSFET  drain  current 
measurements.  The  MOSFET  detection  method  offers  a  number  of  advantages  over  traditional 
piezoresistive  or  capacitive  sensor  elements  because  of  its  small  size,  high  sensitivity,  and  simple 
direct  current  measurement  compared  to  the  complex  piezoresistive  measurements,  as  well  as  its 
full  and  seamless  compatibility  with  direct  monolithic  integration  with  application-specific 
integrated  circuits. 

Moreover,  small  channel  lengths  of  MOSFET  devices  provide  more  localized  stress 
measurements.  MOSFET-embedded  microcantilever  detection  should  allow  for  massively 
parallel  on-chip  signal  sensing,  multiplexing  and  remote-addressability  via  on-chip  integration  of 
RE  (radio  frequency)  elements  as  well  as  photovoltaics  for  local  power  supply  (Figure  28). 


Figure  27B.  Id  versus  Vds  characteristics. 
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Figure  28.  Schematic  illustration  of  potential  of  MOSFET-embedded  microcantilever  array  for  parallel 
and  multiplexed  biomolecular  detection 


4.0  Taking  Lessons  from  T  Lymphocyte  Biology  to  Revolutionize  Detector  Design 


T  lymphocytes  are  the  orchestrators  of  the  adaptive  immune  response  to  pathogen  infection. 

They  detect  the  protein  component  of  antigens  that  have  invaded  cells.  Specialized  cells  called 
antigen  presenting  cells  (APCs)  display  molecular  signatures  of  the  pathogen’s  proteins. 
Specifically,  peptides  (p)  derived  from  pathogenic  proteins  bound  to  protein  products  of  the 
majorhistocompatibility  gene  complex  (MHC)  are  displayed  on  the  surface  of  APCs.  The  T  cell 
receptor  (TCR)  on  the  surface  of  T  cells  can  potentially  bind  these  pMHC  molecules.  The 
surface  of  the  APC  also  displays  abundant  pMHC  molecules  where  the  peptide  was  derived  from 
endogenous  (or  self)  proteins.  TCR  can  also  bind  to  self  pMHC,  albeit  weakly.  Recent 
experiments  carried  out  by  one  of  our  team  members  (Mark  Davis,  Stanford  Univ.)  showed  that 
T  cells  can  detect  as  few  as  ten  stimulatory  pMHC  molecules  in  a  sea  of  30,000  endogenous 
ligands*’.  This  extraordinary  sensitivity  to  antigen  is  not  accompanied  by  frequent  spurious 
activation  events  in  the  absence  of  stimulatory  (agonist)  pMHC  molecules.  Recent  theoretical 
work  carried  out  as  a  collaboration  between  team  members  (Amp  Chakraborty,  UC  Berkeley, 
and  Mark  Davis)  has  shed  light  on  the  molecular  mechanisms  that  enable  sensitive  and  selective 
detection  of  antigen  by  T  cells.  The  goal  of  this  seed  project  was  to  examine  whether  this 
understanding  can  be  harnessed  to  design  pathogen  detectors  that  can  mimic  T  cells  in  their 
ability  to  detect  pathogens.  Toward  this  end,  the  main  goal  of  this  project  was  to  examine 
whether  the  APC  could  be  replaced  by  a  synthetic  surface  and  yet  maintain  the  sensitivity  and 
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selectivity  noted  above.  Specifically,  pMHC  molecules  presented  on  supported  lipid  bilayers 
were  examined  to  determine  whether  or  not  they  can  be  detected  sensitively  and  selectively  by  T 
cells.  Some  theoretical  work  exploring  the  signaling  pathways  of  T  cells  was  also  pursued.  The 
accomplishments  for  this  aspect  of  the  project  can  be  summarized  as  follows: 


1 .  Mixtures  of  pMHC  molecules  presented  by  supported  lipid  bilayers  were  detected  sensitively 
as  measured  by  a  gross  assay  (immunological  synapse  formation). 

2.  Single  molecule  imaging  methods  were  used  to  demonstrate  the  same  as  above  at  the  single 
molecule  level. 


Enhancer  only 


3.  Computer  simulations  were  conducted  to  examine  how  crowding  agents  inside  the  T  cell  can 
influence  signaling,  with  a  view  toward 
exploring  how  signaling  components  put  in  a  Agonist 

synthetic  vesicle  might  work  in  the  absence 
of  crowding  agents. 


Below,  each  of  these  accomplishments  are 
briefly  described. 

Small  amounts  of  stimulatory  pMHC 
molecules  presented  on  synthetic  surfaces 
stimulate  synapse  formation 

Mixtures  of  agonist  and  endogenous  pMHC 
molecules  on  a  supported  lipid  bilayer  were 
presented  using  methods  developed  by  one  of 
our  team  members  (Michael  Dustin,  NYU 
School  of  Medicine)  that  have  been  described 
in  detail  previously'^.  T  cells  were  then 
placed  on  top  of  these  supported  lipid 
bilayers,  and  using  video  microscopy, 
determined  whether  or  not  they  were  able  to 
detect  small  amounts  of  agonists.  When  T 
cells  are  robustly  stimulated,  they  form  an 
immunological  synapse’^.  This  is  a 
recognition  motif  that  consists  of  a  spatially 
organized  collection  of  membrane  proteins. 
The  TCR  and  pMHC  molecules  form  a  tight 
central  cluster  surrounded  by  a  peripheral  ring 
of  adhesion  molecules. 


Enhancer  +  trace  agonist 


Figure  29  shows  that  these  experiments 
demonstrated  that  when  a  supported  lipid 
bilayer  is  loaded  with  as  little  as  0.001  % 
agonist  pMHC  (the  rest  being  a  null  pMHC 


Figure  29.  The  synapse,  a  marker  of  robust 
stimulation,  forms  when  a  large  amount  of  agonist 
is  used  and  for  a  mixture  of  a  small  amount  of 
agonist  with  null  pMHC  (enhancer).  The  enhancer, 
by  itself,  does  not  stimulate  synapse  formation. 
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mimic  of  endogenous  pMHC),  the  T  cell  forms  a  robust  immunological  synapse.  This  indicates 
that,  at  least  in  gross  assays,  T  cells  are  stimulated  by  minute  amounts  of  agonists  in  a  sea  of 
non-stimulatory  pMHC  molecules. 

Single  molecule  detection  of  pMHC  molecules  on  supported  lipid  bilayers  and  assays  of 
calcium  flux 


The  assay  above  only  provides  a  gross  estimate  of  sensitivity.  Therefore,  a  single  molecule 
imaging  methods  developed  in  the  laboratory  of  our  team  member  (Mark  Davis)  which  has  been 
described  in  detail  previously  was  used  to  examine  T  cell  signaling  at  short  times  upon 
interaction  with  supported  lipid  bilayers  containing  mixtures  of  agonist  and  endogenous  ligands. 
An  early  report  of  T  cell  signaling  (calcium  flux)  was  studied  because  this  indicates  whether  or 
not  detection  would  be  fast. 

As  shown  in  Figure  30,  preliminary  results  suggest  that  approximately  3  molecules  of  an  agonist 
pMHC  molecule  presented  on  a  supported  lipid  bilayer  can  be  detected  rapidly  as  determined  by 
fluorescence  imaging  of  calcium  flux. 
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Figure  30.  Single  molecule  imaging  shows  that  roughly  3  molecules  of  an  agonist  pMHC  displayed  on  a 
supported  lipid  bilayer  along  with  a  “sea”  of  endogenous  ligands  can  be  detected. 

The  effects  of  molecular  crowding  on  the  T  cell  signaling  pathway 


As  mentioned  above,  the  ultimate  goal  of  a  full  project  would  include  creating  a  detector  where 
even  the  T  cell  is  replaced  by  a  vesicle  containing  the  minimal  number  of  signaling  components. 
Toward  this  end,  we  wanted  to  explore  how  the  many  inert  macromolecules  in  the  T  cell 
influence  signaling.  We  employed  a  Monte-Carlo  algorithm  to  simulate  a  simplified  version  of 
the  T  cell  signaling  pathway.  We  studied  situations  where  signaling  is  fast  compared  to  the  time 
scale  in  which  the  crowding  agents  move  as  well  as  situations  where  the  inverse  is  true. 

The  main  findings  can  be  summarized  as  follows.  If  signal  transduction  occurs  on  timescales  that 
are  slow  compared  to  the  motility  of  the  molecules  and  organelles  that  constitute  the  crowding 
elements,  the  effects  of  crowding  are  qualitatively  the  same  as  in  a  homogeneous  3 -dimensional 
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medium'^.  In  contrast,  if  signal  transduction  occurs  on  a  timescale  that  is  much  faster  than  the 
time  over  whieh  the  crowding  elements  move,  the  effects  of  varying  the  extent  of  erowding  are 
very  different  when  reactions  occur  both  in  2  and  3-dimensional  space.  For  fast  signaling, 
crowding  agents  attenuate  signaling  and  never  enhances  signaling.  In  contrast,  slow  signaling 
cascades  can  be  both  enhaneed  and  attenuated  by  erowding  agents.  These  results  were  reported 
in  detail  in  a  publieation  in  the  J.  Physical  Chemistry^^.  The  basic  message  from  these  studies  is 
that  the  absence  of  crowding  elements  will  not  be  strongly  influenee  signaling  and  sensitivity 
characteristic  of  T  cells. 

5.0  Summary 

This  effort’s  aceomplishments  can  be  summarized  as  follows:  1 .)  Novel  computational  tools  that 
can  guide  the  design  of  biomolecular  nanomeehanical  systems  were  developed.  2.)  These  tools, 
in  close  synergy  with  experiments,  have  been  applied  to  develop  prototype  devices  using  both 
optieal  and  eleetronie  deteetion  that  are  being  considered  for  commereialization.  3.)  The 
possibility  that  lessons  from  T  eell  biology  eould  be  harnessed  to  create  novel  pathogen  detectors 
were  explored  setting  the  stage  for  future  researeh  opportunity. 
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